Search for the Rare Decay K L — > n°e + e 
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The KTeV/E799 experiment at Fermilab has searched for the rare kaon decay Kl — * 7r°e + e _ . 
This mode is expected to have a significant CP violating component. The measurement of its 
branching ratio could support the Standard Model or could indicate the existence of new physics. 
This Letter reports new results from the 1999-2000 data set. One event is observed with an expected 
background at 0.99 ± 0.35 events. We set a limit on the branching ratio of 3.5 xlO -10 at the 90% 
confidence level. Combining with the previous result based on the dataset taken in 1997 yields the 
final KTeV result: BR{K L -> Ve+e" ) < 2.8 xl0~ 10 at 90% C.L. 

PACS numbers: 13.20.Eb, 11.30.Er, 14.40.Aq 



The decay Kl — ► 7r°e + e _ has long been studied in 
the context of Standard Model CP violation (CPV) and 
has more recently been of interest in certain new physics 
scenarios. 

In the Standard Model, there are direct and indirect 
CPV contributions to the amplitude, plus an interfer- 
ence term [H HI Q The indirect component is known 
from the measurement Q of BR^s — * 7r°e + e~) and 
appears to dominate. The direct component has been 
estimated to be about 3 to 6 xl0~ 12 , and the two CPV 
contributions together give BR(Kl — > n a e + e~)cpv m 
the range 8 to 45 xl0~ 12 . There is also a CP con- 
serving amplitude through 7r°7*7* states which can be 
determined from measurements of Kl — > 7r°77 0,0. 
In recent work, Buchalla, D'Ambrosio, and Isidori |jj 
argue that the CP-conserving contribution is negligi- 
ble. They predict a Standard Model branching ratio 
BR{K L -> 7r°e+e-) ~ 3 x KT U , dominated by CPV, 
with a 40% contribution from direct CPV, through the 
interference term. 

Observation of Kl — * 7r°e + e~ at rates substantially 



higher than Standard Model expectations would signal 
new physics. In a large class of SUSY models, a branch- 
ing ratio enhancement of up to five times the Standard 
Model expectation is considered likely [8(, but values as 
high as 10 -10 are not entirely ruled out. The existing 
experimental limit 9'] has been used to constrain squark 
masses ^3 an d SUSY contributions to the charge 
asymmetry in K^ — > i^i^iT. The implications of a 
specific model with extra dimensions for Kl — * ir e + e~ 
and related processes have been investigated in |T^ . 

The existing experimental upper limit on BR(_/Cl - * 
7r°e+e-) of 5.1 xlCT 10 at the 90% confidence level (CL) 
is based on the 1997 KTeV dataset. In this Letter we 
present an improved limit based on data collected during 
1999-2000. 

At KTeV, 800 GeV/c protons from the Tevatron were 
directed onto a BeO target to create two parallel Kl 
beams. The beams entered a 65m long vacuum tank, 
which defines the fiducial volume for accepted decays. 
Charged particles were detected by two pairs of drift 
chambers separated by an analysis magnet providing a 
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transverse momentum kick of 0.150 GeV/c. Photon ve- 
toes positioned around the vacuum decay region and the 
spectrometer vetoed particles escaping the drift cham- 
bers. The KTeV detector is further described in (ljj. 

Powerful discrimination against charged pions, which 
could fake electrons, was provided by a set of transition 
radiation detectors (TRDs) behind the drift chambers. 
Each of the eight planes was composed of a polypropy- 
lene felt radiator paired with a double-plane multiwire 
proportional chamber containing an 80%-20% admixture 
of Xenon and CO2. TRD cuts resulted in a pion rejec- 
tion factor of about 50:1, as measured in a sample of 
Kl — > n ± e T v decays. These cuts were over 94% effi- 
cient for electrons. A more detailed description of the 
TRD may be found in jlH ]. 

Downstream of the TRDs were the trigger hodoscopes. 
The trigger required hits in the hodoscope planes and 
the spectrometer consistent with the passage of two op- 
positely charged particles. The trigger hodoscopes were 
followed by the Csl electromagnetic calorimeter [lif . 
which had an energy resolution <j(E)/E — 0.45% © 
2%/ ' \J 'E{GeV). Electrons were identified by requiring 
the ratio of the energy measured in the calorimeter (E) 
to the momentum as measured in the spectrometer (p) 
to be consistent with one; this cut rejected about 99.5% 
of charged pions. 

A detailed package of Monte Carlo simulation routines 
was used to study detector geometry and performance, as 
well as various trigger and analysis selection criteria. The 
programs were also used to simulate background events 
and tailor cuts to optimize the signal to background ratio. 

The Kl — > 7r°e + e~ final state consists of two pho- 
tons, which come from the tt° decay, and two electrons. 
Kl — > 7r°e + e~ candidates exhibit the following signa- 
ture: two tracks of opposite charge originating from a 
common vertex, and depositing all of their energy in the 
calorimeter; and two other clusters in the calorimeter, 
which, when taken as photons originating from the ver- 
tex, have a mass consistent with the ir° mass. 

Kl — * k tt d events, where ir^ indicates the pion 
Dalitz decay tt° — > e + e~7, are used to measure the 
Kl flux and normalize the acceptance calculation. This 
mode has a signature similar to Kl — > 7r°e + e~ , with the 
addition of a photon. 

Recorded Kl — > 7r°e + e~ and Kl — > 7T ° 7T d even ts 
satisfied the following trigger requirements. There must 
have been at least two separate track candidates in each 
drift chamber plane. There must not have been hadronic 
showers in the calorimeter, and the event must have 
deposited little energy in the photon vetoes. There 
must have been a minimum number of clusters in the 
calorimeter with energy greater than 1 GeV, as deter- 
mined by the hardware cluster counting system [lfij . For 
Kl — > 7r°e + e~ , this number was four clusters and for 
Kl — > it Tip it was five. 

In the offline event reconstruction and analysis, events 



are required to satisfy further selection criteria. The 
charged tracks must point to calorimeter clusters. To 
identify these tracks as electrons, the ratio of the energy 
of the matched cluster as measured in the Csl (E) to 
track momentum as measured by the drift chambers (p) 
must lie in the range 0.95 < E/p < 1.05. The track posi- 
tions must have sufficient clearance from the Csl edges. 
The decay vertex (Z vtx ) has to be within the vacuum 
decay volume: 96 m < Z vtx < 158 m. The reconducted 
kaon momentum is required to be between 20.3 and 216 
GeV/c. Tracks are required to be well separated (greater 
than 1 cm apart at the first drift chamber) and the open- 
ing angle between the two tracks has to be larger than 
2.25 mrad in the lab frame. 

Further selection cuts for the Kl — > 7r ° 7r D sample 
included requirements on the invariant masses of the 
e + e~7, 77, and e + e~777 combinations, and on the mo- 
mentum transverse to the Kl flight direction, p±. A 
well-reconstructed kaon should have a p± close to zero. 
Using the calculated acceptance and known branching 
ratio for Kl — > t^tvjj decays, the total number of Kl 
decays in the data sample is (349.0 ± 2.8 s t a t ± 21.6 sys t 
± 11.8 B r) x 10 9 . 

Several backgrounds with the e + e _ 77 final state exist 
and can mimic the Kl — > 7r°e + e _ signal. The first source 
of background is Kl —* 7r + 7r~7r° where both charged pi- 
ons shower in the calorimeter and appear to be electrons. 
To remove this background, the mass of the event, under 
the hypothesis that the tracks were pions, is required to 
exceed 520 MeV/c 2 . 

The second source of background is Kl — ► 7r 7r° and 
Kl — » 7r 7r°7r° with one or two Dalitz decays of a 7r° and 
with one or more photons undetected. To ensure that all 
Kl decay products are observed, p\ < 1000 (MeV/c) 2 
is required. Additional background events of this type 
are removed by requiring that the invariant mass of the 
two electrons, m ee , exceed 140 MeV/c 2 . However, there 
are some backgrounds involving two ir^ decays in which 
only one electron and one positron are reconstructed with 
a high mass. These events might also include coincident 
accidental activity. These background events are rejected 
by requiring that m ee be less than 362.7 MeV/c 2 . 

The third source of background is Kl — * 7r ± e^^ where 
the pion fakes an electron by showering in the calorime- 
ter, and photons are radiated by the electron or are ac- 
cidentals. This background is rejected by examining the 
response of the TRDs for both tracks. 

After these cuts are applied, the single largest re- 
maining background is the radiative Dalitz decay Kl — * 
e + e~77 with invariant mass of the two photons, m 77 , 
consistent with the ir° mass. These events come from 
both internal and external bremsstrahlung; both contri- 
butions were studied in |17| . 

In Figure 1, m 77 is plotted against the invariant mass 
of the four-particle system, m ee77 . The m 77 is de- 
termined under the assumption that the photons came 
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FIG. 1: m 77 (charged vertex) vs. m ee77 (neutral vertex) for 
the data after all cuts have been applied except for the phase 
space cuts. The regions appearing in the figure are discussed 
in the text, and signal events in the center box have not been 
plotted. Masses are in GeV/c 2 . 



from the charged vertex, while m ee77 is calculated us- 
ing the "neutral vertex" , found by applying the 7r° mass 
constraint to the photon energies and positions in the 
calorimeter. Better m ee77 resolution is achieved for the 
signal Monte Carlo using the neutral vertex, but this pro- 
cedure gives the incorrect mass for the Kl — > e + e~77 
background, resulting in the diagonal swath in Figure 1. 

There are several distinctive regions in the m 77 vs. 
Tn ee -y~f plane. In order to minimize human bias in the 
determination of the selection criteria, a blind analysis 
was performed. The box was the region covered up until 
cuts were finalized, and spans 130 < m 77 < 140 MeV/c 2 
and 485 < m ee77 < 510 MeV/c 2 . The ellipse in the box is 
the signal region, which spans ^2tr in the Kl — > Tr°e + e~ 
signal Monte Carlo m eeil and m 77 distributions. In the 



m ee77 direction, the ellipse is ± 5.02 MeV/c 2 wide, and 
in the m 77 direction it is ± 2.32 MeV/c 2 wide. The 
rectangular "strip" to the left of the box is dominated 
by backgrounds from Kl — > tt^tt^ and Kl — > ir°irpTTp 
decays with accidental tt s. Missing particles in these 
decays cause the reconstructed to be low. 

Because these backgrounds accumulated in the strip, this 
region is not considered in the background estimation 
described below. 

In the background estimation, the data in Figure 1, 
outside the strip and box regions, is fit to the sum of 
planar parts and the Kl — > e + e~77 sample: 



where g(m ee77 ,to 77 ) is the Kl — *■ e + e~77 distribution 
in the m 77 vs. m ee77 plane. The parameters A t were 
the parameters from the fit. The non-i^ — > e + e~"/"/ 
background was well-modeled with first-order terms. The 
estimated background in the signal ellipse is 38.11 ± 1.67 
events, with 0.27 ± 0.03 event contribution from non 
Kl — » ee77 backgrounds. 

In order to reduce this background, phase space 
cuts are applied to the data. The location of 
these cuts is optimized by minimizing the expected 90% 
C.L. branching ratio limit of Kl — > ir a e + e~ , using the 
Feldman and Cousins methodology. The expected 
branching ratio limit is computed by randomly generat- 
ing a large ensemble of virtual experiments in which the 
known background sources are the only contributions to 
the observed number of events. The effect of the cuts on 
signal efficiency is also accounted for. 

The phase space variables with the best discrimination 
against Kl — > e + e~77 background are | y 1 | and 9 m i n . 
The variable y 1 is the cosine of the angle between the 7r° 
decay axis and the sum of the momenta of the two elec- 
trons, calculated in the center of mass of the photon pair. 
In the signal mode, | y 7 | is nearly uniformly distributed 
because the pion has spin zero, but in Kl — > e + e~77 , 
the distribution is peaked at one. The variable 9 m i n is 
the minimum angle between any photon and any electron 
in the kaon rest frame. It provides good separation be- 
cause in Kl — > e + e~77 , a radiated photon typically has 
a small angle with respect to the electron from which it 
originated, while in Kl — > ir°e + e~ , 9 m in is nearly flat. 
Distributions for | y 1 | and 9 m i n in Kl — * ir°e + e~ Monte 
Carlo and Kl — > e + e~77 data and Monte Carlo appear 
in Figure 2. 

The optimized phase space cut values are 9 m i n > 
0.362±0.017and | y 7 |< 0.745±0.002. These cuts reduce 
the expected background from 38.11 ± 1.67 events to 0.99 
± 0.35 with a signal loss of 27%. The signal acceptance, 
assuming uniform three-body phase space, is (2.749 ± 
0.013)%, giving a single event sensitivity of 1.04 xl0~ 10 . 

When the box in Fig. 1 was opened (Fig. 3), one event 
was observed in the signal ellipse. Taking the background 
level into account, we determine BK(Kl — > 7r°e + e~ ) < 
3.50 xl0~ 10 . Combining this with the previous result 
yields the final KTeV result: BR(K L -> 7r°e+e- ) < 2.8 
xlO- 10 at 90% C.L. 

If instead of a uniform three-body phase space distribu- 
tion for the signal mode, we assume a vector interaction 
model for the direct CPV part of the decay and allow for 
form factors as in 0, we find for the combined 1997 and 
1999 data samples an upper limit of BR(Kl — ► 7r°e + e~ ) 
< 3.4 xl0~ 10 . If the decay K L — * 7r°e + e~ is satu- 
rated by the direct CPV component, we constrain the 
Wolfenstein CKM parameter |??cifM| < 3.3. Although 
other measurements yield a more stringent constraint on 
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FIG. 2: | J/7 | (top) and 8 m i n (bottom) distributions for 
K L -> 7r°e+e~ MC and K L -> e + e"77 data and MC. JO, -> 
e + e~77 events come from inside the swath but outside the 
box. Kl — > 7r°e + e~ MC are from inside the box, and the 
normalization is arbitrary. 
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FIG. 3: m 77 vs. m ee77 in GeV/c 2 for the data after all cuts 
have been applied. The box is open and one event appears 
within the signal ellipse, with a background of 0.99 ± 0.35 
events 



\Vckm\, it is important to make a variety of measure- 
ments in both the kaon system and the B system to de- 
termine if the CKM parameters are consistent. 
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